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ABSTRACT 

We have made detailed calculations of the composition of magnetically driven jets ejected from collapsars, or 
rapidly rotating massive stars, based on long-term magnetohydrodynamic simulations of their core collapse with 
various distributions of magnetic field and angular momentum before collapse. We follow the evolution of the 
abundances of about 4000 nuclides from the collapse phase to the ejection phase and through the jet generation 
phase using a large nuclear reaction network. We find that the r-process successfully operates only in energetic 
jets (> 10^' ergs), such that U and Th are synthesized abundantly, even when the collapsar has a relatively weak 
magnetic field (10'"G) and a moderately rotating core before the collapse. The abundance patterns inside the jets 
are similar to those of the r-elements in the solar system. About 0.01-0.06 Mq neutron-rich, heavy nuclei are 
ejected from a collapsar with energetic jets. The higher energy jets have larger amounts of ^^Ni, varying from 
3.7 xlO"'* to 0.06Mq. Less energetic jets, which eject small amounts of ^''Ni, could induce a gamma-ray burst 
(GRB) a supernova, such as GRB 060505 or GRB 060614. Considerable amounts of r-elements are likely to be 
ejected from GRBs with hypernovae, if both the GRB and hypernova are induced by jets that are driven near the 
black hole. 

Subject headings: Accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances — stars: 
supernovae: general — MHD — methods: numerical — gamma rays: bursts 



1. INTRODUCTION 

During the collapse of a star more massive than 35-40Mo, 
the s tellar core is cons idered to promptly collapse to a black 
hole ('Heger et al."2003'). If the star has sufficiently high angu- 
lar momentum before collapse, an accretion disk forms around 
the hole, and jets has been shown to be launched from the in- 
ner region o f the disk through magnetic processes jProga et al.l 
[2003: Mizun o et aIl2C) 04a b: F ujimoto et al.|2006l) and neutrino 
heating (Nagatak i et al.li2007i) . Gamma-ray bursts (GRBs) are 
expected to be driven by the jets. This scenario for GRBs is 
referred to as the collapsar model (IWooslevlll993h . Assisted by 
accumulating observations that imply an associat ion between 
GRBs and the dea t hs of massive stars (e.g. , Galam a et al.ll 19981 
iHjorth et al.ll2003l:IZeh. Klose. & Hartmann II2004() . this model 
seems to be most promising. 

For accretion rates greater than O.IM0 s"', the accretion disk 
is so dense and hot that nuclear burning is expected to pro- 
ceed efficiently. In fact, the innermost region of the disk, 
thought to be related to GRBs, bec omes neutron-rich through 
electron capture on protons CPruet. Wooslev & Hoffmanll2003l: 
iFujimoto et al.l 12004 l2005ah . Nucleosynthesis inside the out- 
flows from the neutron-enriched disk has been investigated 
with stea dy, one-dime nsional models of the disk and the 
outflows (|Pruet. Thompson & HoffmanI [2004; Fujimo to et al 
2004 l2005bl) . Not only neutron-rich nuclei (.Fujimoto et al 



2004. but also p-nuc lei jPruet. Thompson & HoffmanI |2004 



Fujimoto et al.l l2005bl) has been shown to be produced in 



side the outflows. This co-production of r- and p-elements 
in the outflows has been confirmed with a more elabo- 
rate calcuIa^ion_ofthenucleos in jets from a collap- 
sar dFujimoto et al.ll2007l) . based on two-dimensiona l magneto- 
hydrodynamic (MHD) simulations (iFuiimoto e t al. 2006); the 
calculation, however, was performed only for a collapsar with 
a strong magnetic field (lO'^G) and rapidly rotating core, with 
angular velocity of lOrads"^ In the present study, we recalcu- 



late the chemical composition of jets fro m collapsars, which ar e 
shown to eject magnetically driven jets (iFujimoto et al.ll2006l) . 
to investigate the dependence of jet composition on the collap- 
sar rotation rate and magnetic field. 

Moreov er, the recent discovery of nearby GRB s (060505 and 
060614) ( Fvnbo et al.|[2006l: iGehrels et al.ll2006h that show no 
association with a supernova (SN) has led to the suggestion that 
these GRBs were not produced by a collapsar dGal-Yam et al.l 
I2OO6I: iDeUa Valle etali l2006t IZhang et al.l I2007Ik The ab- 
sence of a SN reveals that very small a mounts of ^^Ni (< 
O.OOlMo) were ejected from these GRBs jFvnbo et al.ll2006l: 
Gehrels et al. 2006), in contrast to the cases of GRBs with SNe; 
the SNe associated with GRB 980425 and GRB 030329, SN 
1998bw and SN 200 3dh, respective l y, ern i tted large amounts 
of 5^Ni (> 0.3Mq) (llwamotoet al.lll998l; iMaeda & Nomotol 
120031) . We show here that the masses of ^^Ni ejected from col- 
lapsars are different from one to another, and as such, the diver- 
sity in the amount of ^^Ni is unlikely to contradict the collapsar 
model. 

In §2, we briefly describe the numerical code used for the 
MHD calculation of the collapsars, or rapidly rotating massive 
stars, the initial conditions of the stars before core collapse, and 
the properties of the jets that are emitted. In §3, we present 
the Lagrangian evolution of the ejecta through the jets and the 
large nuclear reaction network used to follow the abundance 
evolution of the jets. It is also shown that the r-process op- 
erates inside the jets, even when the magnetic field is 10'°G 
and the core is moderately rotating, with an angular velocity of 
2.5rads"'. The masses of ^^Ni and r-elements ejected from the 
collapsars are shown as a function of the energy of the ejecta. 
We discuss the effects of neutrino interactions and the energy, 
liberated through nuclear reactions on nucleosynthesis in the 
jets and the dynamics of the collapsars in §4. Finally, we sum- 
marize our results in §5. 
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2. MHD CALCULATIONS OF COLLAPSARS 

We have carried out two-dimensional, Newtonian MHD cal- 
culations of a collapse of rotating, massive (40Mq) star whose 
core is assumed to promptly collapse to a black hole. Here 
we present our numerical models and the results of the MHD 
simulations, in particular, the production and properties of 
the jets, which are important for nucleosynthesis. For de- 
tails of the MHD mo dels and results, the reader is referred to 
iFuiimoto etal] (120061) . 

2.1. Input Physics and Numerical Code 

The numerical code for the MHD calculations employed in 
this p aper is based on the ZEUS-2D code ( Stone & Norman 
1 19921) and is the same as used in iFuiimoto^t alj (l2007l) . We 
have extended the code to include a realistic equation of state 
(EOS) (Kotake et al. 2004) based on relativistic mean field the- 
ory (iShen et al.l 11998). For the lower density regime (p < 
lO^gcm"^), where no data are available in the Shen EOS ta- 
ble, we use another EOS (Blin nikov et al. l fT996l) . We consider 
neutrino cooling processes. The total neutrino cooling rate is 
evaluated with a simplified neut rino transfer model based on 
the tw o-stream approximation ( Di Matteo. Perna. & Naravaiil 
l2002l) . with which we can treat the optically thin and thick 
regimes of the neutrino reactions approximately. We ignore 
resist ive hea ting, the properties of which are highly uncer- 
tain (iProga et al. 2003). We note that viscous heating is not 
taken into accounts. We assume the fluid is axisymmetric and 
has mirror symmetry about the equatorial plane. Spherical co- 
ordinates, {r,9,(f)) are used in our simulations and the com- 
putational domain extends over 50 km < r < 10,000 km and 
< 6* < 7r/2 and is covered with 200(r) x 24(6*) meshes. Fluid 
is freely absorbed through the inner boundary at 50km, which 
mimics the surface of the black hole. The mass of the black 
hole is continuously increased by the infall of gas through the 
inner boundary. We mimic strong gravity around the black hole 
in ter ms of a pseudo-Newtonian potential dPaczvriskv & Wiitai 
I1980I) . 

2.2. Initial Conditions 

We set the initial profiles of t he density, t emper ature, and 
electron fraction to those of the iHashimotol ( Il995h spherical 
model of a 4OM0 massive star before the collapse. The radial 
and azimuthal velocities are set to zero initially and increase 
as a result of the collapse induced by the central hole and self- 
gravity of the star. The computational domain extends from the 
Fe core to an inner O-rich layer and encompasses about 4-Mq of 
the star The boundaries of the Si-rich layer between the Fe core 
and the O-rich layer are located at about 1 800 km ( 1 .88M0) and 
3900km (2.4M0), respectively. We adopt an analytical form for 
the angular velocity 17 of the star before the collapse: 



as in previous studies of coUapsars dMizuno et al ]l2004illbl) and 
SNe (iKotake et al.ll2004) . Here flo and Rq are parameters of our 
model. We consider three sets of (fio, ^0): (lOrads"', 1000km) 
(rapidly rotating core) (2.5rads"', 2000km) (moderately rotat- 
ing core) and (0.5rads"', 5000km) (slowly rotating core). We 
note that for these values of ^Iq and Rq, the maximum specific 
angular momentum is about 10^^ cm^ s"\ which is comparable 
to that of the Keplerian motion at 50km around a 3M0 black 
hole. Therefore, the centrifugal force can be larger than the 



gravitational force of the central black hole, and the formation 
of a disk like structure is expected near the hole. The initial 
magnetic field is assumed to be uniform and parallel to the ro- 
tational axis of the star We consider cases with initial mag- 
netic fields (Bo) of 10^ 10>°, and lO'^G. It should be noted that 
the magnetic pressure is much smaller than the gas pressure 
initially, even for Bq = lO'^G. We have performed MHD sim- 
ulations for nine models, labeled by R8, RIO, R12, M8, MIO, 
M12, S8, SIO, and S 12, where the first character, R (rapidly ro- 
tating core), M (moderately rotating core), or S (slowly rotating 
core), indicates the set of flo and Rq and the numeral, 8, 10, or 
12, equals logB()(G). Note that the MHD simulations for col- 
lapsars with a moderately rotating core (M8, MIO , and Ml 2) 
are in addition to those presented in lFuiimoto et al.l d2006) . 

2.3. Properties of Jets 

We briefly describe th e results of our MHD calcula- 
tions (Fujimot o et alj|2006l) . The Lagrangian evolution of the 
physical quantities inside the jets, which is important for nucle- 
osynthesis, is described in §3.1 in detail. 

We find that jets are magnetically launched from the cen- 
tral region of the star along the rotational axis; after mate- 
rial reaches the black hole with a high angular momentum of 
^ lO'^cm^ s"', a disk is formed inside a surface of weak shock, 
which appears near the hole and propagates outward slowly, be- 
cause of the centrifugal force. The magnetic fields are chiefly 
amplified by the wrapping of the field inside the disk, so that the 
toroidal component increases and dominates over the poloidal 
component. Eventually, the jets can be driven by the tangled 
magnetic field lines at the polar region near the hole. 

We find that jets are launched for six models, RIO, MIO, 
SIO, R12, M12, and S12. Their properties are summarized in 
Table 1 . Columns (2) through (4) list the initial conditions for 
the models. Column (5) gives the time at the end of each run, 
t = tf, after the onset of collapse, or f = 0. Column (6) gives 
fjet, the time when the jets pass through 1000km. Columns (7), 
(8), (9), and, (10) represent the mass and magnetic, kinetic, and 
internal energies of the jets, respectively. We note that the jets 
are more energetic for the coUapsars with a moderately rotating 
core (MIO and Ml 2) than for R12. Larger amounts of gravita- 
tional energy are liberated in the central parts of coUapsars MIO 
and M12 than in R12. This is because the centrifugal force of 
the core for R12 is larger than those for MIO and M12. The ac- 
cretion rates onto the black hole through the accretion disk are 
larger for MIO and M12 than for R12; the rates are 0. IM0 s"' 
for MIO and M12 and - O.OIM© s"' for R12 at the epoch of the 
launching of the jets. 

3. NUCLEOSYNTHESIS IN JETS 

We next examine nucleosynthesis in jets from coUapsars, 
models RIO, MIO, SIO, R 12, Ml 2, and SI 2, based on the re- 
sults of the MHD simulations ('Fujimot o et al.ll2006l) . We first 
describe the Lagrangian evolution of the ejecta through the jets 
in detail and then proceed to nucleosynthesis inside the jets. 

3.1. Lagrangian Evolution of Ejecta through Jets 

In order to calculate the chemical composition of the jets, 
we need the Lagrangian evolution of physical quantities, such 
as the density, temperature, and velocity of the jet m aterial. 
We adopt a tracer-par ticle method dNagataki et al.lll997) . as in 
iFujimoto et al.l (l2007h . to calculate the Lagrangian evolution 
from the Eulerian evolution obtained from our MHD calcula- 
tions. Particles are initially placed between the Fe core and an 
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Fig. 1 . — Maximum densities, pmax, and temperatures, Tmax, of jet particles for models with Bq = lO'" G (left) and Bq = lO'- G (right). The circles, triangles, and 
squares indicate models with rapidly, moderately, and slowly rotating cores, respectively. 



inner O-rich layer. The total numbers of the particles A^^ are 
set to be 50,000, 1000, 2000, 1000, 2000, and 5000 for RIO, 
MIO, SIO, R12, M12, and S12, respectively, with which we 
can follow the ejecta through the jets appropriately. We show 
how different values of Np change the masses and abundances 
of the ejecta through the jets in 33.31 The number of particles 
in a layer are weighted according to the layer's mass. We find 
that 21, 105, 18, 59, 214, and 113 particles are ejected through 
the jets for models RIO, MIO, SIO, R12, M12, and S 12, respec- 
tively. Hereafter, we refer to these as jet particles. 

We continued the MHD calculations until t = tj (Table 1). 
After tf, we assume that the particles are adiabatic and expand 
spherically and freely. Therefore, the velocity, position, den- 
sity, and temperature of each particle are set to be v(f) = vo, 
r{t) = rf + v(t)(t-tf), pit) = pf{rf/r{t))\ T(t) = Tf(rf/r(t)), re- 
spectively, where vq is constant in time and set to Vf. Here Vf, 
rf, pf, and 7/ are the velocity, position, density, and tempera- 
tu re of the particle at tf , respectively. It should be noted that 
in lFujimoto et alj ( l2007l) we discussed how the assumption re- 
garding the expansion affects the abundance changes in the ex- 
pansion phase of the ejecta for a low particle in model R12, 
and we concluded that the particle composition depends only 
weakly on this assumption. Here, is the electron fraction at 
9 X 10'' K. 

For the strongly magnetized collapsars (R12, M12, and S 12), 
the jet particles are initially lo cated inside the Fe co re or just 
above it, as shown in Figure 1 of lFuiimoto et al.l ( l2007h for R 12, 
while the particles are initially placed in the O- or Si-rich lay- 
ers for the mildly magnetized collapsars (RIO, MIO, and SIO). 
As the particles fall toward the black hole, their densities and 
temperatures rise, and then they decrease when t he particles are 
ejected through the jets, as shown in Figure 2 of lFuiimoto et al.l 
(l2007h for R12. Figure [T] shows the maximum densities and 
temperatures of the jet particles, pmax and Tmax- Particles with 
relatively low density (< lO'gcm"^) have p^ax a nd r^ax sim- 
ilar t o those of ejecta from Type II SNe ( Thi elemann et al.l 
119981) . On the other hand, for MIO, R12, and M12, some 
jet particles have extremely high densities and temperatures 
(Pmax > lO^gcm"^ and Tmax > 10'"K), which cannot be real- 
ized in the ejecta from a Type II SN. These particles are dense 
and hot enough for electrons to be captured on protons, so that 
the particles become neutron-rich, as we show below (Figure|2]). 
Hence, the r-process is expected to operate in the jets for MIO 
and M12, as well as R12, in whic h the r-process has al ready 
been shown to operate successfully dFuiimoto et al.ll2007b . 

3.2. Nuclear Reaction Network 



As shown in Figure [T] some ejecta can attain to tempera- 
tures higher than 9 x lO^K near the black hole, at which point 
the material is in nuclear statistical equilibrium (NSE). The 
abundances of material in NSE have simple analytical expres- 
sions (Clayton 196 8), which are specified b y p, T, and the elec- 
tron fraction, as in lFuiimoto et al] ( l2005a[) . The electron frac- 
tion changes through electron and positron captures on nuclei. 
It should be emphasized that the changes in the electron frac- 
tion through neutrino interactions can be ignored in the ejecta, 
as we explain in 34.11 

In the relatively cool regime of T < 9 x lO^K, NSE breaks 
down and the chemical composition is calculated with a large 
nuclear reaction network (network B in Nishimura et al. 200g). 
The network includes about 4000 nuclei from neutron and pro- 
ton up t o fermium, whose atom ic number Z = 100 (see Ta- 
ble 1 of iNishimura et alj l2006l) . The network contains re- 
actions such as two and three body reactions, various de- 
cay channels, and electron-positron capture (for details, see 
iFuiimoto et al.l 120071) . Experimentally determined masses and 
reaction rates are adopted if available. Otherwise, theoreti- 
cal nuclear data, such as nuclear masses, rates of two body 
reactions, and /3-decays, are calculated with a mass formula 
based on the extended Thomas-Fermi plus Strutinsky integral 
(ETFSI) (iGorielv. Tondeur & Pearsonll200ll) . Spontaneous and 
/3-delayed fission processes are taken into accounts in the net- 
work. We note that the dependence of abun dances of r-elem ents 
on the mass formula are discussed in lFuiimoto et al.l ( l2007l) for 
model R12. 

3.3. Composition of Collapsar Jets 

Once we obtain the density and temperature evolution of a 
jet particle, we can follow its abundance evolution during in- 
fall and ejection, through a post-processing calculation using 
the NSE code and the nuclear reaction network described in the 
previous subsection. We integrate the abundances of all jet par- 
ticles weighted by their masses to obtain the composition of the 
jets from the collapsars. 

The maximum densities and temperatures of the jet parti- 
cles pinax and Tmax (Figure [1]), are good indicators of the nu- 
clear yields of the ejecta, which are experienced explosi ve nu- 
clear burning (see Fig. 8 of iThielemann et alj 1199 8*). For 
Tmax < 5 X lO^K, the particles are abundant in a-elements, such 
as ^'^Si, ^^S, ^''Ar, and "^''Ca, due to incomplete Si burning or ex- 
plosive O burning. On the other hand, for Tmax > 5 x lO^K the 
particles undergo complete Si burning, which leads to the syn- 
thesis of a large fraction of ^''Ni. The jet compositions for mod- 
els RIO and SIO are therefore similar to those of ejecta from 
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the Fe- and Si-rich layers of Type II SNe. It should be empha- 
sized that the ejected masses for a-elements lighter than Fe are 
underestimated because ejecta from an outer O-rich layer (> 
10,000 km) are not taken into account. 

Figure |2] shows Yg with respect to p,^.^^ for the jet particles. 
Here is the electron fraction at 9 x 10^ K of the the parti- 
cles. We note that all ejecta have Y^ ^ 0.5 for RIO and SIO 
and Ye > 0.467 for S12, while there exist many particles with 
lower values, Y, < 0.4, for MIO and M12, as well as R12. The 
r-process operates in these low-Yg or neutron-rich, jet particles. 
Figure [3] shows the time evolution of the electron fraction, en- 
tropy per baryon, dynamical timescale, which is evaluated as 
r/vr, and temperature for the most neutron-rich jet particle in 
MIO. Here Vr is the radial velocity of the particle. As the par- 
ticle falls in toward the black hole, it passes through a shock 
surface 0.44 s), and the temperature and entropy are en- 
hanced as a result of shock heating. Despite the short dynam- 
ical timescale (< 10~^s) during the infall near the black hole, 
the electron fraction decreases through rapid electron capture 
on protons due to the high density and temperature. This rapid 
capture leads to fast neutrino cooling, which lowers the entropy 
efficiently 0.48 -0.68 s). After the ejection of the particle 
through the jet, which may enhance the entropy of the parti- 
cle ('^ 0.7 -0.9 s), the temperature drops to less than 5 x 10^ K 
and neutron capture on nuclei proceeds efficiently because of 
the relatively slow ejection (fdyn > 0.1 s). The entropy of par- 
ticles with low Ye is low [~ (5- 10)^6, where is the Boltz- 
mann constant, and the dynamical timescale is relatively long 
0.01 -0.3 s), as shown in FigurelH in which the entropy per 
baryon and dynamical timescale are shown for T = 9 x 10^ K. 
The neutron-rich particles (7^ < 0.2) are likely to have a high 
neutron-to-seed ratio (> 100) and a low num ber of fission cy- 
cles, as discussed bv lMever & Brown 1 (fT997l see their Figs. 5 
and 6), although they considered cases with much higher en- 
tropy [-(50-500)/tB]. 

Figure|5]shows the mass distributio n with respect to Y^ in the 
ejecta for MIO and M12 (see Fig. 9 of lFuiimoto et alj|2007l for 
R12). We find that the masses of particles with lower Yg (< 0.4) 
are 0.031 and O.OISM© for MIO and M12, respectively, which 
are comparable to the value for R12 (0.018Mq). Note that the 
ejected masses, M^, are 0.222 and O.129M0 through the jets for 
MIO and M12, respectively (Table 1). Therefore, large amounts 
of r-elements are ejecte d through the coUapsa r jets for MIO and 
M12, as in model R12 (Fujimoto et al. 2007). 

The compositions of jets as a function of mass number. A, are 
shown in Figure|6]for MIO (left) and M12 (right). The jets from 
the collapsars of MIO and M12 have abundance profil es simi- 
lar to that of solar r-elements, just as for R12 (Fuiimo To et al.l 
l2007l) . although there are details in the profiles that are rather 
different from that of the solar r-elements. We note that ap- 
preciable amounts of U a nd Th are produced in the jets for 
MIO and M12, as in R12 (Fujimoto et al]|2007h . We also note 
that the collapsar jets have abundant r-elements compared with 
Fe, as we shall see clearly in Figure [8] below. In the collap- 
sar jets, [Eu/Fe] and [Ba/Fe] are ~ 4 and ~ 3, respectively, 
much larger than the values for the metal poor star, HE2148- 
1247, which shows large enhancements of r-elements ([Eu/Fe] 
= 2.0 and [Ba/Fe] = 2.4) ( Cohe n et al, 2003). Here [A/B] is 
log{YA/Ys)-log(YA/YB)Q for elements A and B with number 
fractions Ya and Yg- 

For MIO and M12, the collapsar jets have large amounts of 
light p-nuclei (^'^Se, ^^Kr, ^"^Sr, and ^^Mo) and heavy p-nuclei 



("■'in, '^^Sn, and '^^La). These p-nuclei are much more abun- 
dant than is found in core-collapse S Ne and are produced w ith- 
out i-process seeds, as in model R12 dFujimoto et al.ll2007l) . 

The ejected mass slightly increases for larger A^,,. The values 
of Mej are 0.222, 0.121, and 0.084 Mq for MIO, M12, and R12, 
respectively, with Np = 1000. When we use 2000 tracer par- 
ticles, these increase to 0.245, 0.129, and 0.124 Mq for MIO, 
Ml 2, and R12, respectively. The mass distribution with respect 
to Ye also depends weakly on A'^ for Y^ < 0.4, and thus the abun- 
dance profiles of r-elements change slightly for different values 
of Np. Figure|7]shows the mass distribution with respect to Yg in 
the ejecta of M12 with half the number of particles (A',, = 1000). 
Compared with A^^ = 2000 case (Figure |5] right), ejecta with 

~ 0.3, in which nu clei are synthesized m ainly with A around 
130 (see Figure 4 of dFujimoto et alj|20()7l) ) are deficient. The 
r-elements with A ~ 130 are less abundant than for the case 
with Np = 2000 (Figure |6] right). We conclude that the abun- 
dance profiles of the jets do not depend strongly on A^p. We 
note that the profiles also depend weakly on the nuclear mass 
formula and th e numerical resolution of the MHD simulations, 
as discussed in lFujimoto et al.l d2007h . 

3.4. Masses and Energies of Collapsar Jets 

The jets are abundant in ^^Ni because of the complete Si 
burning, as described in ^3.31 (Figure [TTi. Figure [8] shows the 
mass of ^''Ni with respect to the energy of the jets, where the 
energies are the sum of the kinetic energy (£'j;)ej and the inter- 
nal energy (£',)ej of all the jet particles. The plus sign, open 
square, open circle, open triangle, cross, and asterisk indicate 
^''Ni for models RIO, R12, MIO, M12, SIO, and S12, respec- 
tively. We find that the ^''Ni masses are roughly proportional 
to the energies of the jets. High energy jets (MIO, R12, M12, 
and S12) could be observed as GRBs with normal SNe if the 
collapsars are located near the Earth. On the other hand, low 
energy jets (RIO and SIO) might be observed as GRBs with- 
out SNe , like the recently observed nearby GR Bs, 060505 and 
060614 dFvnbo et alJl2006l:lGehrels et alj|2006h . even if the jets 
are accompanied by GRBs located near Earth. 

The ejected masses of ^^Ni in our simulations are 
lower than those in the hydrodynamic simulations of 
iMacFadven & WooslevI dl999 ), in which a large amount of 
^^Ni is expected to be ejected through winds driven by vis- 
cous heating. In contrast, in our MHD simulations, ^^Ni 
is ejected through the jets, not through winds, which disap- 
pear in the simulations. The viscous heating is implemented 
with the a-prescription in MacFadyen & Woosley's hydrody- 
namic code while such heating is not taken into account in 
our MHD simulations. As a result of viscous heating, the en- 
tropy of the winds is greater than (20- 30)A;b per baryon and is 
much higher than that of the jets in our MHD simulation (Fig. 
[3]). We note that win ds have not appeare d in MHD simula- 
tions by other groups dProga et alj 2003; N agataki et al.ll2067l) 
or in the hy drodynamic simulation with a low a ^ 0.001 by 
IMacFadven & Wooslev ( 1999). The winds are therefore gener- 
ated by an "a-viscosity" and may be artificial. 

It should be emphasized that the velocities of the jets are less 
than 0.2c, where c is the speed of light, in our MHD simula- 
tions. The jets are therefore too slow to produce GRBs. We 
note, however, that the low-density polar region near the black 
hole in the collapsars is a possible candidate for the produc- 
tion site of GRB jets. This is because the region is located near 
the inner region of the accretion disk, in which a large mag- 
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Fig. 2.— Electron fraction, 7^, at 9 X IO'K vs. maximum density of jet particles, pmax, for models with Bq = lO'^G (left) and Bq = lO'-G (right). Circles, 
triangles, and squares indicate jet particles for models with rapidly, moderately, and slowly rotating cores, respectively. 
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Fig . 3 . — Time evolution of the electron fraction, entropy per baryon, dynamical timescale (r/vr), and temperature for the most neutron-rich jet particle in model 
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Fig. 4. — Dynamical timescale (squares) and entropy per baryon (circles) vs. Ye for the neutron-iich jet particles in model MIO at 9 X lO' K. 



netic energy > 10^' ergs is sustained and throu gh which mass 
accretion takes place at rates above 0.05Mq s"^ ( Fuiimoto et al.l 
l2006h . The energy liberated in the disk can be greater than 9 x 
10^' ergss"' even for the low efficiency (=0.1) of energy libera- 
tion. In addition, the energy of neutrinos, which emana tes from 
the disk at rates higher than 5 x 10^' ergss"' dFuiimoto et al.l 
12006) . is partly transfered to the polar region from the disk, and 
acts to heat the polar region dNagataki et al.ll2007h . 

Moreover, collapsar jets may have two components: ultra 
relat ivistic, baryon-poor j ets and non relativistic, baryon-rich 
jets dSoderberg et alj2006l) . SNe associated with GRBs are per- 
haps induced through an aspherical, jet-like expl osion, which 
produces non relativistic, baryon-rich jets (.Maeda & Nomotol 



l2003h . The SNe accompanying GRBs suggest that the jets have 
two components. The jets in the present study correspond to 
the non-relativistic, baryon-rich jets. 

The masses of r-elements are also revealed in Figure [8] with 
the filled square, filled circle, and filled triangle standing for 
models R12, MIO, and M12, respectively. The r-process is 
found to be operate only in energetic jets (> 10^' ergs ), which 
have considerable amounts of ^''Ni. The energetic explosion 
leads to the ejection of not only a large amount of ^^Ni but 
also considerable amounts of neutron-rich material, on which 
the r-process operates, near the black hole. If jets are realized 
with higher energies in a collapsar, they probably have larger 
amounts of ^^Ni as well as appreciable amounts of r-elements. 
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Fig. 7. — Same as Figure|5]but for model M12 with half the number of particles (Np = 1000). 



The jets can be observed as GRB associated with a hypernova, 
in which abundant r-elements are possibly produced. There- 
fore, GRBs with hypernovae, such as GRB 980425 and GRB 
030329, could eject appreciable amounts of r-elements in addi- 
tion to a large amount of ^^Ni, if the hypernova is induced by 
jets that are driven near the black hole. 

We terminated the simulations at f = f/ (Table 1), which is 
much shorter than the typical duration ~ 10 s of a long GRB. If 
we perform the simulations for a much longer time, the ejected 
masses and energies will increase. Figure |9] shows time evo- 
lution of mass ejection rates, Mej, and integrated mass, Mgj, 
ejected through a spherical surface with a radius of 1000 km, 
after the launch of the jets for models MIO (0.7-1.3 s) and 
M12 (0.21-0.32 s). Mass ejection takes place at rates of 0.2- 
0.5 Mq s"' for MIO, while the jets for M12 have much higher 



ejection rates, up to 5M0S . The average ejection rates from 
fjet to tf (Table 1) are 0.42 and 1.17 MqS"! for MIO and M12, 
respectively. If these high ejection phases were to continue for 
10s, the ejected masses could be 4.2 and 11.7 Mq for MIO 
and Ml 2, respectively. It should be emphasized that the GRB 
duration indicates the continuous ejection of ultra-relativistic, 
baryon-poor jets, not baryon-rich ones. The mass fraction of 
^^Ni in the ejecta is about 0.3 for MIO and M12. The masses 
of ^''Ni, therefore, can be comparable to those of hypernovae 
for the ejection of 1-2M0 material through the collapsar jets, 
which is realized in the ejection of the jets during a period of 
2.5-5 s and 1-2 s for the high ejection phases in MIO and Ml 2, 
respectively. These ^^Ni-rich ejecta could have large amounts 
of r-elements, ^ 0.1 -0.2Mo, because of the large abundances 
(^ 0.1 in mass) of the elements for MIO and M12. We need. 
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however, much longer term simulations to determine the dura- 
tion of the ejection of baryon-rich jets and the masses of ^''Ni 
and r-elements. 

4. DISCUSSION 

4. 1 . Effects of Neutrino Captures 

We first discuss the effects of neutrin o capture on the compo - 
sition of the coUapsar jets. As shown in Fuii moto et alj (l2007h . 
jet particles with lower have more abundant heavy nuclei. If 
the neutrino capture on nucleons is taken into account, the elec- 
tron fraction of the particles changes compared with that in the 
present study, so that heavy nuclei may be less abundant. Tak- 
ing into account the change in electron fraction due to neutrino 
captures on nucleons, we follow the evolution of the electron 
fraction of the jet particles for model MIO, in which r-elements 
are abundantly ejected through the jets (Figure |6] left). Then 
we estimate the difference in between cases with and without 
neutrino capture for all the jet particles. We assume that neu- 
trino emission is isotropic for simplicity, although it actually is 
aspherical because neu trinos are chiefly emitted from the dense 
and hot disk (Fig. 9 of 'Fujimoto et al.' 2006). We evaluate the 
capture rates of electron neutrinos and anti-neutrinos, Aj^^ and 



Ap,, as 

A,, =4.83F„L,,,5i£.7'-i (2) 
X^^=4.S3YpU^,5iE.-Jrl (3) 

where Y„, Yp, Lp,,5i, E^^, Ec;^, and r^ are respectively 

the number fractions of neutrons and protons, the luminosities 
of electron neutrinos and anti-neutrinos in units of 10^' ergs~\ 
the average energies of electron neutrinos and anti-neutrinos 
in MeV, and the distance between a particle and the center of 
the coUapsar in units of 10^ cm (Oian & Woosley 1996). We 
set Ei,^ and /i^^ to 10 and 15, respectively. We have evalu- 
ated the neutrino luminosity with the tw o-stream approxima- 
tion dPi Matteo. Perna. & Naravanll2002l) in our MHD simula- 
tions. However, the luminosity is the sum of neutrino luminosi- 
ties over all neutrino flavors, and thus L,y^,5i and Lp^.si are not 
provided by the two-stream approximation. We t herefore esti- 
mate L^^ 51 and according to Appendix B of lRuffert et al.l 
( 1996), with the density and temperature profiles obtained from 
the MHD simulations. We note that the evaluated total lumi- 
nosities are comparable with the above two methods. We find 
that .51 and .51 are less than 3.6 and 3.1, respectively, for 
MIO, and that Y^ increases compared with the case without the 
neutrino capture, especially for particles with Ye < 0.2. How- 
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ever, the differences in Yg are small, up to 0.03. There exist 
particles with Yg < 0.15, and the r-process operates in the jets 
to produce U and Th abundantly. We note that for some low Y^ 
particles, Y^ decreases compared with the case without neutrino 
capture. 

Moreover, we followed the evolution of the electron fractions 
of jet particles for model MIO, with L^^.si and L,7^,5i set 5 times 
larger than the luminosities evaluated with the results of the 
MHD simulation. This is because when the inner boundary of 
the computational domain for the MHD simulations, rin, is set 
to 10 km instead of 50 km as in the pres ent study, t he neutrino 
luminosity known to increase slightly dFuiimoto e t al. 2006). 
Aspherical neutrino emission from the disk may increase the 
neutrino intensity on the jet particles effectively. We find that 
the difference in Y^ is as much as 0.12 and that there are sev- 
eral particles with low Y^ < 0.2, in which U and Th are abun- 
dantly synthesized. The total mass of the particles amounts to 
0.01 IMq. We conclude that the r-process probably operates in 
the collapsar jets to synthesize considerable amounts of U and 
Th even if we take into account neutrino capture on nucleons 
and aspherical neutrino emission. 

In addition to the changes in electron fraction due to neutrino 
capture, neutrino interactions are importa nt for the dynamics o f 
collapsar jets through neutrino heating ( Nagataki et alj|2007l) . 
The heating from neutrino captures and annihilations between 
neutrinos and anti-neutrinos may increase the energy of the jets. 
The amount of ^^Ni ejected through the jets is likely to increase 
and could produce hypernovae with appreciable amounts of r- 
elements. 

4.2. Ejecta from the Outer Layers of Collapsars 

In the present study, we have calculated the masses and com- 
position of the jets ejected from the inner region (< 10,000 km) 
of the collapsars. The ejecta from the outer layers (> 10,000 
km) have a large amount of material and are abundant in ele- 
ments lighter than Fe. 

It should be emphasized that comparable amounts of ^''Ni 
may be ejected through the jets from outer O-rich layers (> 
10,000 km), which are not covered by the computational do- 
main of the present MHD simulations. We have evaluated 
M(^^Ni)out for ejecta from the core-collapse SN of a 40 M© star, 
which is the initial model for our MHD simu lations, using the 
spherical explosion model jHashimotoll 1 995h . Here M(^^Ni)out 
is the mass of ^''Ni ejected from the outer layers (> 10,000 km). 
We set the explosion energy of the SN to be 1.0 x 10^' ergss"', 
which is comparable to the energy of the ejecta from the col- 
lapsar of model MIO. We find that 0.1 15M0 of ^*'Niis ejected 
from the outer layers. If we take into account an aspherical jet- 
like explosion, M(^''Ni)out is estimated to be ^ 0.04Mq. This is 
because, for MIO, the jets have a relatively large opening angle 
of ~ 30 degree. We note that M(^^Ni)out is comparable to the 
mass of ^*Ni from the inner layers < 10,0 00 km. 

As discussed in Fujimo to et al.l (l2007h . the /7-process in O- 
rich layers is expected to operate in colla psars, as in the p - 
process layers (PPLs) of core-collapse SNe dRavet et"al]|1995h . 
This is because i-process seeds for the p-process exist in the 
layers, as a result of the weak ^-process operating during He 
core burning of the stars, and the jets propagate in the layers 
to heat them to high enough temperatures for the p-process to 
operate efficiently. Therefore, the masses of /^-nuclei can be en- 
hanced by means of the /^-process. The /^-nuclei ejected from 
collapsars with high energy jets represent the sum of those in 



the ejecta from oxygen-rich layers, with composition similar to 
that in the PPLs, and those in the jets from the inner core, which 
abundantly contain p-nuclei that are deficient in the SNe, such 
as ''^Mo, ii^^In, "^Sn, and '^^La. 

4.3. Effects of Nuclear Energy Generation 

We finally discuss how the energy that is generated through 
nuclear reactions affects the dynamics of collapsars. We eval- 
uate the ratio, /?nuc, of the specific energy liberated in nuclear 
reactions, enuc, to the specific energy, e, of a particle. Here e 
is obtained from the results of the MHD simulations and enuc is 
calculated as 

e„ue = 9.65 X 10'^ [ Q> ' (4) 

where the sum is performed over all the nuclei in the reaction 
network and Qi and 5Yi are the mass excess (in MeV ) and the 
change in the number fraction of the / th nuclide, respectively. 
We obtain (5F, from the abundance evolution of the particle. We 
find that |/?nuc | is always less than 0.05 in jet particles, other than 
neutron-rich jet particles with Yg < 0.4. Figure [TOl shows the 
time evolution of the magnitude of /?niic for a jet particle with 
^ 0.1 in model R12. Solid and dotted lines represent cases 
with positive and negative Rnuc, respectively. As the jet particle 
falls in near the black hole (f < 0.2 s), heavy elements in the 
particle are chiefly destroyed to "^He and then to protons and 
neutrons by photodisintegrations, which are endothermic. On 
the other hand, during ejection through the jets (0.2 s < f < 1 s), 
nucleons in the particles recombine to "*He and then heavier nu- 
clei, chiefly neutron-rich nuclei, by means of the r-process. The 
nuclear energy liberated by the recombination heats the particle 
because of the exothermic nature of the recombinations. Af- 
ter the freeze-out of the r-process due to neutron consumption, 
abundant neutron-rich unstable nuclei begin to decay through 
sequences of /3-decays (t > 0.4 s). The liberated energy is com- 
parable to the particle energy, but \R„uc\ is always less than 
0.487. We conclude that the nuclear energy is unlikely to affect 
the dynamics of the collapsar jets, significantly, as in supernova 
ejecta with a neutron-to-proton ratio of 10 (see. Fig. 3 of iSatol 
[1971 and in ejecta with an ele ctron fraction of 0.12 from neu- 
tron star mergers (see, Fig.2 of iFreiburghaus et al ] IT999l) . 

5. SUMMARY AND CONCLUSIONS 

We have calculated the composition of magnetically-driven 
jets ejected from collapsars, or, rapidly rotating massive 
(40Mq) stars, based on long-term, magnetohydrodynamic sim- 
ulations of the stars during core collapse. The magnetic fields 
of the collapsars prior to the collapse are set to be uniform and 
parallel to the rotational axis. We consider three sets of angular 
velocity distributions, with rapidly, moderately, and slowly ro- 
tating cores, and three cases with magnetic fields of 10'^, lO'", 
and 10'^ G. We follow the evolution of the abundances of about 
4000 nuclei in the jets from the collapse phase to the ejection 
phase through the jet generation phase with the aid of an NSE 
code and a large nuclear reaction network. We summarize our 
conclusions as follows: 

1 . The collapsars eject appreciable ^^Ni, in amounts com- 
parable to or somewhat less than that from normal core- 
collapse SNe. The ejected ^^Ni masses, which range 
from 3.7 x 10"'* to O.O6M0, through the jets are roughly 
proportional to the jet energy. 
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2. Less energetic jets that eject small amounts of ^''Ni 
could induce GRBs without a supernova, such as GRB 
060505 and GRB 060614. 

3. The r-process operates only in energetic jets (> 
lO^^ergs), even if changes in the electron fraction due 
to neutrino captures are taken into account. U and Th 
are synthesized in jets even for a collapsar with a mod- 
erately rotating core of 2.5rads"' and magnetic field of 
lO'^G. 

4. The abundance patterns of the energetic jets are similar 
to those of the solar r-elements, while low-energy jets 
have compositions similar to the ejecta from core col- 
lapse SNe. 

5. Considerable amounts of r-elements could be ejected 
from GRB with a hypernova, if both GRB and hyper- 
nova are induced by jets that are driven near the black 
hole. 

6. Energy released through nuclear reactions is not signif- 
icant for the dynamics of the collapsar jets. 

In the present study, we have calculated the masses and com- 
position of the jets ejected from the inner regions (< 10,000 
km) of the collapsars. The ejecta from the outer layers (> 
10,000 km) contain a large amount of mass and are abundant in 



elements lighter than Fe. As discussed in 34.21 the amount of 
^^Ni in the ejecta from the outer layers is comparable to that 
from the inner layers, and considerable amounts of p-nuclei 
are ejected through the jets from the outer layers. We need to 
take into account the ejecta from the outer layers to estimate the 
masses and abundances of elements lighter than Fe, as well as 
p-nuclei. This will be our next undertaking. 

Relativistic effects related to black hole spin could induce 
vary energetic jets at a later time when the rotation of the black 
hole becomes rapid as a result of continuous mass accretion 
with high angular momentum through than accretion disk. Am- 
plification of toroidal magnetic fields by the magneto-rotational 
instability might also produce more energetic jets if we per- 
form a three-dimensional MHD simulation, rather than two- 
dimensional simulations as in the present study, in particular, 
for models with low-energy jets, such as RIO and SIO. In 
addition, the standing accretion shock instabilit y, whose im- 
portance has been recognized in SN explosion C Blondin et al.l 
2003) might change the energies and masses of the jets. 
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TABLE 1 
Jet Properties 



model 


Bo 


Oo 


Ro 




fjet 




(^m)ej 






RIO 


lOio 


10 


1000 


2.62 


2.58 


0.0010 


2.89e-4 


0.0274 


0.0840 


MIO 


10'" 


2.5 


2000 


1.30 


0.77 


0.222 


0.787 


2.206 


7.926 


SIO 


10'" 


0.5 


5000 


1.34 


1.30 


0.0053 


0.045 


0.138 


0.606 


R12 


1012 


10 


1000 


0.36 


0.20 


0.083 


0.097 


4.58 


2.79 


M12 


1012 


2.5 


2000 


0.33 


0.22 


0.129 


0.783 


24.5 


9.44 


S12 


1012 


0.5 


5000 


0.28 


0.25 


0.033 


0.014 


3.57 


5.27 



Note. — The model parameters, Bq, ^Iq, and R^, are shown in units of G, rads" , and 
km, respectively. The calculations are stopped at time tf. The jets pass through a point 
1000km from the central black hole at time ?jet. The mass ejected through the jets, Mq, 
is in units of Mq. The magnetic, kinetic, and internal energies carried away by the jets, 
(£'m)ej, iEk)ej, and iEi)ej, are expressed in units of lO^^erg. 



